NACA EM E50D25 


•RESTRICTED 


ft ’ 


Copy 

RM E50D25 


- Min-< 



' 


RESEARCH MEMORANDUM 


ANALYTICAL INVESTIGATION OF FLOW AND HEAT TRANSFER IN 
1 COOLANT PASSAGES OF FREE -CONVECTION LIQUID-COOLED 

1 Q 3 

I i TURBINES 

j 3 

■g cd By E. R. G. Eckert and Thomas W. Jackson 

& , <% K 

i ; j l Lewis Flight Propulsion Laboratory 

j I ‘ Cleveland, Ohio 


.*> 1 "i 


2 1 


CLASSIFIED DOCUMEHT 


ThU document cocUIni elusiOad tnfirmition 
affecting a» Hatlcnai Defers# erf th* United 
gtitu within th* tcawpinc of tkw Espfccog* Act, 
USC 50:31 and St. Its tncsmlnS» or tbw 
raswlstirs at It* eon tant* Ek isf nmunr to an 
um'jUjortMd ptmi la prokltaltad lv In, 
InfcimiUcB. ao etaaalQed ma y b> imparted 
otTj to p»rioB> tn On military ud wnl 
Mnkti a t Eh* Unttsd SHtu, ipfrcprlit* 
olvlllu officers auf cmpLaynis st the Fsdsnl 
Gev»TOE.«* wt» t*v* s lwgBlir.it> lnt*r#sL 
tbwrwi*, md to United States ciUssns at knewrn 
loyalty ird diweretio* wto at necessity must bw 
Inf armed tfcarsof . 


i'- A - A UHRARY 

LK*n0.r'i il^AL LAaawrar 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

July 18, 1950 


UMCLASSiriFf 4, 



KACA EM E5GD25 


T i *n 

^■BSTOvTiai/ 


UNCLASSIFIED 


NATIONAL ADVISORY COMMITTEE FOE AEEONADTICS 
BESEAECH MEMORANDUM 

ANALYTICAL INVESTIGATION OF FLOW AND HEAT TRANSFEE IN COOLANT 
PASSAGES OF FREE -CONVECTION LIQUID -COOLED TURBINES 
By E. R. G. Eckert and Thomas V. Jackson 

SUMMARY 

Gas turbines provided with free -convection cooling have 
rotating blades with cylindrical holes that are closed on top 
of the blade and connected with the fluid supply at the blade 
root. In the design of such blades, an important question 
arises as to what minimum diameter a hole of a certain length 
can be drilled without endangering the. circulation of the cool- 
ing liquid in it. An approximate answer to this problem is 
obtained by calculating the thickness of the heated boundary 
layer that builds up around the cylindrical surface of a hole 
with large diameter. 

Numerical evaluations of the results of these calculations for 
a sample turbine and comparisons with test results show that in 
small holes such as would be required in the leading and trailing 
edges of turbine blades, the circulation and the cooling effect is ap- 
preciably impaired. The circulation and the cooling effect in these 
holes can be considerably improved when a connection is provided 
on top of the blades between the small holes and the larger holes, 
which can be located in the central portion of the turbine blade. 

In the boundary-layer calculations, the effect of the Coriolis 
forces on the movement of the liquid in the rotating holes was 
neglected. The magnitude of the Coriolis forces relative to the 
centrifugal forces was calculated and the influence of the secondary 
movement set up by the Coriolis forces on the cooling effect was 
estimated by comparison with test results on inclined surfaces in 
gravitational free-convection heat transfer. 


INTRODUCTION 

Liquids are more effective coolants than gases because of 
their better conductivities. Gas turbines can therefore he 
operated with higher gas temperatures and, accordingly, with better 
efficiency and greater power per unit gas flow when they are cooled 
with liquids. 
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The most critical parts to he cooled in gas turbines are the 
rotating blades. In using liquids to cool the blades, the very 
large free -convection currents that are generated by the centri- 
fugal forces can be used. This cooling method was first proposed 
by E. Schmidt (reference 1). 'Experiments at the NACA Lewis lab- 
oratory proved its effectiveness. Figure 1 shows this cooling 
arrangement for a turbine blade according to the proposal of 
Schmidt. The cooling liquid fills up the cylindrical holes of 
circular oross section inside the blade. The holes are closed 
on top and connected at the bottom with the fluid supply. 

When the blades are heated by hot gases, the heat penetrates 
through the solid material of the blades to the surface of the 
holes, where it heats a thin layer of liquid adjacent to the 
surface. When the liquid layer is heated, the density in the 
layer becomes smaller than the density in the cooler liquid core 
near the center of the hole. In the rotating turbine, the fluid 
in the hole is subjected to a centrifugal acceleration equal to 

ra?. (All symbols used in this report are defined in appendix A.) 

The centrifugal forces pro? per unit volume are smaller in the 
heated layer (boundary layer) because of its smaller density p 
and oause this layer to flow toward the axis of rotation. The hot 
fluid leaving the hole in this way has to be replaced by cool fluid 
flowing outward in the core . The velocity profile in a oross 
section of the hole must therefore have a shape similar to that 
shown in the central hole of figure 1. 

The circulation of the liquid in a hole will be good as long 
as the oross section, of the core is large as compared with the 
cross section of the heated boundary layer. When the diameter of 
a hole is small as compered with its length, then the space left 
for the core will become too small and the circulation will 
diminish. For holes with very small diameters the flow may break 
down entirely. 

On the other hand, it le necessary to arrange as many holes 
as possible in the blade in order to minimize the thickness of 
solid material through which the heat must flow. If the heat has 
to flow over a long distance within the solid material, the blade 
cannot be oooled effectively. An analytical investigation of the 
problems arising in connection with this cooling method was con- 
ducted at the NACA Lewis laboratory and is presented herein. This 
analy sis investigates : (1) the smallest diameter hole that can be 
made without endangering the circulation of the liquid, and (2) 
methods of improving the circulation in a small -diameter hole. 
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METHOD OF ANALYSIS 

A hole of very large diameter as compared with its length 
is first considered. In this hole, the cross-sectional area 
of the core of the liquid is very large compared with the 
cross-sectional area of the heated layer adjacent to the wall 
and consequently the velocity of the fluid in the core is very 
sma ll . By decreasing the diameter of the hole, the flow in the 
heated layer and its thickness will change; however, as long as 
the cross section of the core is considerably larger than the 
cross section of the layer, the change will he negligible. Only 
when both areas are comparable in magnitude will the circulation 
start to decrease appreciably. Consequently, if it is possible 
to calculate the thickness of the heated layer for a very large 
hole, it is possible to estimate at what diameter the circulation 
will start to decrease appreciably. 

As long as the length of the hole is small compared with 
the diameter of the turbine wheel, so that the radius in the 
expression for the centrifugal forces may be replaced by some 
average value, the conditions for the flow in such a hole are 
practically the same as the conditions of free -oonvect ion flow 
in a vertical stationary cylinder that is filled with liquid and 
whose walls are heated. The bottom of the cylinder is closed and 
the top is open and connected to a liquid supply. A liquid film 
around the inside of the wall will be heated and rise by its 
buoyancy forces . Cold liquid from the supply on top will replace 
the warm liquid film leaving the cylinder. The centrifugal force 

prc£ that causes the circulation of the cooling liquid in the 
turbine blade is replaced by the gravitational force pg (per unit 
volume) in the equipment just described. 

The flow conditions in a hole with a very large diameter sure 
the sane as on a vert ical -plane plate. Considerable information 
exists on gravitational free -convection flow on a vertical wall. 

It is characterized by a dimensionless value, the Grashof number, 
in the same way that forced-convection flow is characterized by 
the Eeynolds number. In comparing the gravitational free- 
convection flow with the centrifugal flow, the acceleration 
due to gravity g has to be replaced by the centrifugal accele- 
ration no 2 in the expression for the Grashof number. When the 
Grashof numbers in both systems have the same value, the results 
obtained theoretically or experimentally on a gravitational flow 
can be transferred to the flow in the turbine -blade holes. 
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A bou n da r y layer "builds up on the surface of the vertical wall 
in free -convection flow in much the same manner as in forced - 
convection flow. The flow within the boundary layer is 1 ami nar or 
turbulent, depending on the magnitude of the Grashof number. The 

transition takes place at a Grashof number between 10® and 10?"®. 

The thickness of the laminar boundary layer, as well as the 
temperature and velocity profiles within this layer, and the heat 
transfer were measured by E. Schmidt and W. Beckmann (reference 2) 
and compared with values that they calculated with the help of 
E. Pohlhausen. The Grashof numbers in the holes of liquid -c ooled 
turbine blades, based on the length of the holes, are in the 

approximate range of lO-^ to lCp-^. The boundary layer is therefore 
expected to be turbulent over almost the whole length of the holes. 
Some measurements exist on the heat transfer in the turbulent 

region, but they cover only a range of Grashof numbers from 10 9 to 

10- 1 - 2 (references 3 and 4) . Nothing, however, is known on the thick- 
ness of the turbulent boundary layer or on the velocity and tempera- 
ture profiles in this layer. 

Insufficient information is available on turbulent forced 
flow to calculate exactly velocity profiles and boundary-layer 
thicknesses. In addition, practically no information is available 
on the mechanism of free -convection flow in the turbulent range. 

It is therefore possible to make only approximate calculations 
for the boundary -layer thickness in this region. The calculations 
using the ’S&rsahx momentum equation of the boundary layer 
(reference 5) and an analogous heat-flow equation are presented in 
the section DETERMINATION OF TURBULENT FREE -CONVECTION BOUNDARY 
LAYER. Simi lar calculations were made by Squire (presented in 
reference 6) for the laminar free -convection boundary layer and 
his results satisfactorily checked the results of the exact theory 
and of experiments. 

So far only the centrifugal forces in the fluid that fills 
the holes have been considered. It Is known, however, that in 
studying the motion relative to a rotating system, two kinds of 
forces have to be considered, the centrifugal and Coriolis forces. 
The Coriolis forces set up a secondary movement in the liquid 
circulating In the holes of a rotating turbine blade that is 
superimposed on the movement generated by the centrifugal forces. 

The relative magnitude of the centrifugal and Coriolis forces is 
considered and the secondary movement described in a subsequent 
section of the paper. 
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DETERMINATION OF TUKBULENT FREE -CONVECTION BOUNDAET USER 
Derivation of B oundary -Layer Equations 

If a stationary, plane, vertical nail is heated to a tempera- 
ture higher than its* surroundings, the layer of fluid adjacent to 
the vail is heated hy conduction from the vail. In this vay, 
"buoyancy forces are generated that cause this layer to flow in an 
upward direction. A boundary layer is huilt up that begins vith 
zero thickness at the lover end of the vertical vail and increases 
in thickness in the upward direction. The boundary layer becomes 
turbulent, depending on the critical Grashof number, at a certain 
distance from the lower end of the vail . The distance measured 
vertically from the lover end of the plate is called x and the 
distance normal to the plate is y. In order to determine the 
boundary -layer thickness for steady state, a small stationary 
volume element in the turbulent region of the boundary layer is 
considered. Figure 2 shows this volume element. The dimensions 
of the element are dr along the plate and l normal to it. 

The length l should be larger than the boundary-layer thickness 
5. Considering two-dimensional flow, the dimension of the volume 
element normal to the plane of figure 2 may be considered to be 
unity. The upward velocity of the fluid in plane 1-1 at a dis- 
tance y from. the surface of the vail is u. Then the mass flov 
through a small area with a width dy is pu dy and the flow of 
momentum in the x -direct ion is pa? dy. The mom entum flow in the 
r -direction entering the volume ele m ent through plane 1-1 is 

f pu 2 dy 
J 0 

In progressing to plane 2-2, the momentum flow changes by 



In general, the mass flow entering the volume element 
through plane 1-1 is different from the mass flow leaving it 
through plane 2-2, which means that fluid enters or leaves the 
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volume element through the plane parallel to the wall at a 
distance Z. Inasmuch as It is assumed that the velocity in the 
x -direction can he neglected outside the boundary layer, no 
momentum in the x -direction is carried through the plane. 

The rate of change of momentum must he in equilibrium with the 
forces acting in the x -direction on the fluid within or on the 
surface of the volume element considered. A shearing stress T w 
acts on the wall. The force connected with this stress is t w dx. 

No shearing stress occurs on the surface of the volume element 
that is parallel to the wall and at a distance I from it because 
outside the boundary layer the velocity in the x-direction is zero. 

According to the boundary -layer theory, the pressure change 
can be neglected along any normal to the surface. A constant pres- 
sure difference dp therefore exists between planes 1-1 and 2-2. 
This pressure difference gives a force on the volume element of 
magnitude Z dp. In addition to the forces on the surfaces, there 
is a force, due to the weight of the fluid within the volume element 

Z 

pg dx dy 


Summing up all the forces and equating them to the change in 
momentum flow gives the momentum equation 



— / (pu 2 dy) dx = Z dp - 

dx 


pg dy dx - t v dx 


Outside of the boundary layer there is no flow and the pres- 
sure difference between planes 1-1 and 2-2 is therefore balanced 
by the weight of the fluid between the planes 


dp = p 5 g dx 

Multiplying both sides of the last equation by Z and ohanging the 
right side of the equation to the integral form gives 


Z dp 



p 5 g dx dy 
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Introducing the preced in g equation into the momentum equation 
gives 


1 


d_ 

dx 



(pu 2 dy) dx 


l 



(p s -p) dx dy - r w dx 


By introducing the expansion coefficient defined by the 
equation 


P v 5 

the first term on the right side of the preceding equation can 
be transformed to 


g dr / pp (t-tg) dy 


By designating the difference between 
distance y and the temperature tg 
layer by &, the following expression 


the temperature t at the 
outside of the boundary 
is obtained: . 


g dx 



Pp0 dy 


In the applications considered, the changes in the density 
and in the expansion coefficient are smalJ compared with their 
mean values. Therefore, both quantities can be assumed constant 
in the preceding expressions. The momentum equation for the free 
convection boundary layer therefore be comes 
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A similar equation is set up for the heat flow through the 
volume element in figure 2 . The heat oarried with the fluid 
through plane 1-1 is 


8 


where the enthalpies Cp0 are measured from the temperature out- 
side of the "boundary layer. The specific heat is considered 
constant . 

The heat carried out of plane 2-2 differs from that carried 
into plane 1-1 "by 


pCpU0 dy 



The difference in the heat flow through planes 1-1 and 2-2 must 
come from the surface and therefore the heat flow leaving the 
plate per unit time and area is 


% 


gpc. 


l 

h/ 

Jo 


vl& dy 


( 2 ) 


Equation (2) is the heat -flow equation for the free -convection 
boundary layer. 


Solution of Boundary-Layer Equations 

According to von K&rm^n (reference 5), equations (1) and (2) 
may he solved by assuming certain shapes for the velocity and 
temperature profiles within the boundary layer and an empirical 
relation for the friction force on the wall. By introducing 
expressions for the velocity and temperature profiles into the 
preceding momentum and heat-flow equations, the integrations in 
the equations can be performed. Two total differential equations 
result, from which two unknown values (for example, the boundary - 
layer thickness and the maximum velocity within the boundary layer) 
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can 'be obtained as functions of x. Naturally, the results of 
such a calculation are "better if the assumed temperature and vel- 
ocity profiles closely approximate the actual temperature and 
velocity profiles. It was found, however, that the method is 
quite insensitive in respect to the shapes assumed. The difficulty 
in calculating the turbulent free-convection boundary layer in 
this way is that no information whatsoever exists on the shape of 
either profile. It is only known that the approximate shapes must 
be similar to those shown in figure 3. The exoess temperature 6 
has the value © w at the wall and decreases steadily towards the 
value zero, which is reached at a distance 5 from the wall. The 
velocity u has the value zero at the wall and outside of the 
boundary layer and therefore increases to a maximum value 

with increasing distance from the wall and then decreses again to 
zero. 


The velocity profile in the forced -convection turbulent 
boundary layer is approximated quite well by the expression 



(3) 


The shape of the velocity profile in the free-convection boundary 
layer mast be different because the velocity is zero outside of 
the boundary layer as well as on the wall. In the neighborhood 
of the wall both velocity profiles will probably be similar. 
Equation (3) is therefore modified by a factor that brings the 
velocity to zero at the outer border of the boundary layer 


u • u i (|) 1/7 ( x - 1 ) 

When the maximum velocity is determined by differentiation, 

it is found that 
% 


U]_ * 1.54 i i^t (5) 

Equation (4) represents the velocity profile that is used in 

the first calculation. The profile is shown as ( — — ) in fig- 

VW/i 

ure 4 together with two other velocity profiles . 
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In his calculation of the turbulent boundary layer on a rotat- 
ing disk, von KaiWtn used the same velocity profile and found good 
agreement of the calculated friction values with experimental 
results (reference 5). The calculation of the turbulent free- 
convection boundary layer is repeated with the other velocity pro- 
files in figure 4 in order to study their effects. 

The temperature profile in the turbulent forced-convection 
flow for Prandtl numbers not too far from 1 is s imi lar to the 
corresponding velocity profile and it will probably be similar in 
shape for free -convection flow. The following equation is there- 
fore used for the temperature profile: 



The shearing stress t w on the wall in turbulent forced- 
convection flow is given in reference 5 as 

l/4 

T w = 0.0225 pu 2 ^iL\ (7) 

By assuming that the same law holds for turbulent free -convection 
flow in the immediate neighborhood of the wall and considering 
that the velocity profile (equation (4)) is identical with that 
represented by equation (3) for very small distances from the wall, 
equation (7) can be transformed to 


T 


w 


0.0225 pu-j^ 



By using equation (4) for the velocity distribution, the 
integral of the first term in the momentum equation (equation 1) is 




0.207 u^S 


0.491 u 5 
max 


In the preceding equation, the upper limit of the integral was 
changed from l to 2 = 1 because equation (4) holds only for 
y < S, whereas for y > 5 the velocity is zero. 
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The integral of the Becond term, in equation (1) is 

X 9aj ' v i 

The integral in the heat -flow equation (equation (2)) is 

^ . u l 6 o w r (I) 1 / 7 (1- I)[i-(|) 1/7 ] d(|) 

Jq 


= 0.0681 u-^S 


= 0 .105 u ^, T e v s 


Substituting the values of the preceding integrals in the 
momentum and heat-flow equations and considering that the heat- 
transfer coefficient H can he substituted for the specific 
heat flow q w = H0 V gives 

eE ('W 2 ®) = I ef» v « - 0 - 0479 %ai 2 C®> 


H ' fe = °- los 8P0 p fe (^ 8 ) 


(9) 


Solution of equations (8) and (9); Method I . - In order to 
solve equations (8) and (9), additional information on the heat- 
transfer coefficient H is necessary. Jakob (reference 4) 
derived from a compilation of experiments the empirical relation 

1/3 


Nu = 0.129 (Gr Pr) 


( 10 ) 
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Beoause the Nusselt number Is proportional to z and the Grashof 
number Is proportional to x 3 , the foregoing equation gives a 
heat -transfer coefficient that is independent of z. With this 
relation, the differential equation (9) can easily he integrated. 
Considering the houndary condition that for x * 0, 5*0, the 
result is 


u max 6 


Hz 

0.105 gpop 


(H) 


The differential quotient in equation (8) can. he transformed hy 
differentiation as follows: 

df (w* 8 ) ■ u max a! [v) + u max 5 

Introducing the values from equations (9) and (11) for the product 
Ujuj^S and its differential quotient in equation (8) gives the 
following differential equation for the velocity u ^y i 


Au max 

dx 


0.254 


%az 


^tax 

x 


0.00582 


( SPOp \ 5 ^ 4 Umax 2 * 1 / 4c 

\ H ) ^574“ 


( 12 ) 


No analytical solution was found for this equation, which was 
therefore solved graphically, as described in appendix B. The 
results can he approximated with good accuraoy hy the equation 


Be 


max 



0.377(Gr Pr) 0 * 500 


Pr 


-0.65 


(13) 


Up to now, the houndary-layer thickness has heen represented 
hy the value 5, which is the distance from the well where the 
velocity (as approximated hy equation (4)) and the temperature 
(according to equation (6)) reaoh zero values. A more character- 
istic measure for the houndary-layer thickness in general use is 
the so-callepl displacement thickness, which has to he defined for 
free -convection flow hy the equation 


8 
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which. is the thickness the boundary layer would have if the -vel- 
ocity within it were constant and equal to the value Umi T and 
the volume flow were the same as in reality. 

By using the velocity ratio as determined from equations (4) 
and (5), the displacement thickness "becomes 


6* = 0.628 8 (14) 

In equation (11) the heat -transfer coefficient H and the 
maximum velocity u^-r can "be replaced "by the Grashof and Prandtl 

numbers "by means of equations (10) and (13). The houndary-layer 
thickness 8 can he ohanged to the displacement thickness 8* hy 
means of equation (14), which results in 

^ = 2.04 (Gr Pr)“°‘ 167 Pr“ 0 * 35 (IS) 


Solution of equations (8) and (9)t Method 2 . - Because the 
graphical solution is a comparatively tedious process, another 
approach to solve the differential equations (8) and (9) is tried. 
It will he investigated how closely the equations can he satisfied 
with the assumption that u^^ and 8 are proportional to some 


power of the distance x. 

If 


* 

Umax « C u a“ 

(16) 

and 

8 » CqX? 1 

(17) 


are introduced, equations (8) and (9) "become 


0.491 C u 2 C 5 (2mfn)x 2m+n “ 1 * | gpe^Qi 11 - 



(7/4)m-(l/4)n 


(18) 


% 


= H 


and 


0.105 gpc p C u C 5 (m+n) 


(19) 
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Inasmuch, as these equations must he valid for any value of x, all 
the exponents of x in any of the equations must he equal. From 
the momentum equation, values of m and n are obtained: 

m ** l/2 

n = 7/10 

It can he seen from the heat -flow equation that the heat- 

transfer coefficient must he proportional, to x^-/® in order to 
fulfill this equation. 

By -writing 


H « CjjX 1 / 5 


( 20 ) 


and introducing this expression into the heat-flow equation. 


0.126 gpCpCyCg « Cg 


( 21 ) 


is obtained. This equation can he solved for 


u 0.126 gpCpCg 


( 22 ) 


When the value of C u from equation (22) is introduced into 

the momentum equation, keeping in mind that m = 1/2 and n « 7/10, 
the momentum equation becomes 


I sew - ^ - 0 


SPOp j 


(23) 


which is a cubic equation for Cg that characterizes the houndary- 
layer thickness . 
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Xt is usually assumed that in free -conveot ion flow the Husselt 
number is a function of the product of the Grashof and Prandtl num- 
bers. In order to fulfill this condition and equation (20), the 
Nusselt number for turbulent free -convection flow must be represented 
by 

Nu » C(Gr Pr) 2 / 5 (24) 

The range of G-rashof numbers for which test results are available 

(10^® < Gr < 10 12 ) is too small to decide whether equation (10) or 
(24) is the better one. The test points that are presented in 
figure 25-3 of reference 4 and figure 129 of reference 7 agree 
even slightly better with the slope of equation (24). (See fig. 5.) 
Comparison of this equation with the preceding test results gives 
for the constant C the numerical value 0.0210. Therefore, 

Ifu - 0.0210 (Gr Pr) 2 / 5 (25) 

If the displacement thiokness 8* is made dimensionless by 
dividing by a length x and is then introduced into equation (23) 
together with the dimensionless parameters 


Ir-fit 

k 

Hu = — “ 0.0210 (Gr Pr) 2 / 5 
k 

that characterize the flow and heat transfer in free convection, 
the following equation is obtained: 

- 0.0734 (Gr Pr)" 0 * 2 Pr" 1 * 0 ^ - 


0.00415 (Gr Pr)" 0 * 3 Pr" 0,75 * 0 


(26) 
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Tbs local heat-transfer coefficient is needed for the calculations 
presented herein, whereas the test points in fig. 5 are heat- 
transfer coefficients averaged, over the plate length. A slight 
difference exists between both values ; however, inasmuch as the 
Grashof numbers for the test points are too low to determine the 
exact numerical value of the constant in equation (24), the con- 
stant 0.0210 will be used. Because test results are usually 
plotted aB Nu against (Gr Pr), the produot (Gr Pr) was 
introduced into equation (26) . 


Equation (26) was numerically solved for several values of 
(Gr Pr) . The results are presented in figures 6 and 7 as 


( 8 , 


Inasmuch as in both logarithmic charts the values for the dimension- 
less boundary -layer thickness are practically straight lines, the 
following expression approximates the result of the cubic equation 

in the indicated range 10^ < (Gr Pr) < 10^ very closely: 


» 0.296 (Gr Pr)" 0 * 100 Pr -0 * 411 (27) 

The second unknown value was the maximum velocity within the 
boundary layer. This value is represented by the Reynolds number 

From equations (16) and (22), is determined as a func- 

tion of Cg. Replacing Cg by means of equations (20) and (25) 
gives 

*W,I ■ °- 105 -TTT® 

Pr ' ' 

By introducing equation (27) for the boundary -layer thiokness, 
Ee max,I = °* 355 ( Gr Pr) 0,500 Pr“°‘ 589 (28) 


is obtained 
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*■ In. figure 8 the values Ke max and S*/x, as determined by 

hoth methods of solving the differential equations, are compared. 
It can he seen that the differences are not excessive in the 
entire range of Grashof numbers Investigated, when the lack of 
information on turbulent free -convection flow is considered; in 
no case do they exceed 25 percent. 


Solutions of Boundary-Layer Equations Using Other 
Telocity Profiles 

The Barman method used in the preceding section proved com- 
paratively insensitive to the assumed shape of the velocity and 
temperature profiles in forced -convection flow. In order to check 
this insensitivity for the free -convection flow, the calculations 
were repeated using two other expressions for the velocity profiles. 

“ n ; - u i (|) 1/7 [ x - ( i )] 2 ( 29 > 


*1 - 1 * 69 Umax' 



U 1 = 2 * S0 u max 


The shapes of these velocity profiles, together with the first 
assumption and the temperature profile, are shown in figure 4. The 
calculations resulted in the following equations : 

“ 0.0315 (Gr Pr)“°* 2 Pr -1 *° - 


0.00218 (Gr Pr)“°* 3 Pr -0,75 = 0 


(31) 
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£i\ 5 - 0.0308 (Gr Pr)" 0,2 Pr" 1,0 

, Z /IH 


0.00583 (Gr p r )“°* 3 Pr “ 0,75 =0 (32) 

For the dimensionless displacement thickness of the "boundary 
layer, the equations can "be approximated in the range of Grashof 

numbers from 10 ^-® to 10^ 3 "by 

= 0.205 (Gr pj.) -0.100 pj.-0.358 ( 33 ) 


(§1) * 0.235 (Gr Pr )” 0 * 100 Pr -0 - 328 (34) 

\ X /III 

The Reynolds numbers characterizing the .maximum velocities 
can be approximated by 


Ee max,II 3 0.418 (Gr Pr ) 0 * 500 Pr -0 * 642 (35) 

Be max,UI - 0.384 (Gr Pr ) 0 * 500 Pr " 0 * 672 (36) 


Discussion of Results 

The boundary-layer thicknesses and the maximum Reynolds num- 
bers, as- determined with the three velocity profiles, differ only 
in the constant and slightly in the exponent of the Prandtl number. 
The average value of the constant for the boundary -layer thickness, 
based on the three values that were obtained with the different 
velocity profiles, is 0.245 and the three individual constants 
deviate by ±20 percent from this average value. The corresponding 
average for the constant in Ee max 1 b 0.386 and the individual 
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values deviate by ±8 percent from this average value. So it can 
he seen that the influence of the shape of the assumed profiles 
is not large. Of course, there may he an additional deviation of 
the calculated values from, the real ones, which is caused hy the 
uncertainty of the value for the shearing stresses. The calcula- 
tions can he improved in this direction as soon as more information 
on this value in turbulent free -convection flow is available. 

It is quite interesting to see hy which factors the coef- 
ficients m and n in equations (16) and (17) are influenced. 

A check of the calculations leading to these equations shows that 
the value m = l/2 follows from the assumption that the houndary- 
layer thickness and the maximum velocity Increase with same power 
of x. The same value for the coefficient m is valid in the 
laminar range. The value for the coefficient n also depends on 
the law for the shearing stress T w . a value of n may he deter- 
mined hy assuming that the heat-transfer coefficient is independent 
of x. 

The calculation in this manner gives the result that the shear- 
ing stress T v has to depend on the power -l/2 of the Keynolds 

number uy/v in order to fulfill this demand. This large deviation 
frcm the law for the shearing stress in forced -convection turbulent 
flow (equation (7)), where the corresponding power is -l/4, seems 
improbable. The shape of the velocity or temperature profile does 
not influence m or n. 

Another assumption involved In the calculations should he 
mentioned. In assuming that the boundary-layer thickness S is 
proportional to s? 1 , It is implied that the turbulent boundary 
layer begins with the thickness zero at the lower end of the 
plate. In reality, first a l aminar boundary layer occurs and then 
the turbulent boundary layer develops within a certain distance x 
out of the laminar boundary layer. The same assumption that Is made 
here is ordinarily used for calculating the turbulent boundary 
layer for forced flow along a plate and proves quite successful 
(reference 8) . In addition, the relative range of the laminar 
boundary layer at the high Gras hof numbers encountered In turbine 
operation is so small that an error In the correct starting point of 
the turbulent bo undar y layer influences the result of the calculations 
to a very small degree. 
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Comparison with Experiments 

A one -stage gas turbine with liquid free -convection cooling 
was tested at the HACA Lewis laboratory (reference 9) . The results 
of these tests may be used to make a first check on the calculations . 
The tested turbine wheel had a mean radius r_ of 5.72 inches. 

SI 

Each blade of the wheel had five holes with a length L of 2 .5 
inches. Three holes were 0.125 inch in diameter, one was 0.090, and 
one was 0.060 inch in diameter. The circumferential velocity of 
the blades at the mean radius was approximately 700 feet per second. 
This velocity gives 1,050,000 feet per second per second for the 
centrifugal acceleration. This high acceleration, which is 32,000 
times the acceleration due to gravity, explains the large Grashof 
numbers attained in these holes. The cooling water had a tempera- 
ture of approximately 110° F. The blade -wall temperature was 425° P. 
By using the preceding values and a film temperature of-200° F for 
the property values, the Grashof number is found to be 1.25 X 1CP-4. 
Considering the Prandtl number to be 1.74, the Grashof number to be 

1.25 X 1CP-^, and the kinematic viscosity to be 3.16 x 10“® square 
feet per second, introducing these values into equations (27) and 
(28), and solving for 5* and u mar gives 6* » 0.0214 inch and 

u mai — 58 feet per second. If the diameter of a cooling hole is 

determined in such a way that the water velocity has the same value 
in the core of the fluid as in the heated boundary layer, the value 
for this diameter is found to be 0.146 inch. In a hole of such a 
diameter, the heat -transfer coefficient has surely already decreased 
as compared with a big-diameter hole. A comparison of this value 
with the diameters of the holes that were used in the test turbine 
leads to the conclusion that the circulation of the cooling water in 
these holes must have been quite restricted. 

In the tests, the average heat-transfer coefficient for all 
holeB was measured. The corresponding Nusselt numbers are plotted 
in figure 5. The same figure shows the Nusselt numbers calculated 
with equation (24) and with equation (10) given by Jakob. Both 
Nusselt numbers refer to the heat -transfer coefficient in a hole 
with a large diameter. Seme test results are also shown in figure 5. 
They were taken from figure 25-3 of reference 4 and figure 129 of 
reference 7. It can be seen that equation (24) and Jakob's equation 
fit the. test results equally well and that the Nusselt values given 
by both equations are very much larger than the Nusselt numbers 
measured in the test turbine; so the conclusion, which iB drawn from 
the calculations , that the holes in the test turbine were too small 
in diameter for proper oiroulation, is confirmed. 
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Inasmuch as it is impossible to avoid holes with small dia- 
meters in turbine blades, the question arises as to how the cir- 
culation can be improved in these holes. This problem will be 
dealt with in the following section. 

A brief calculation reveals that vaporization will never 
occur in a hole if its diameter is large enough for good cir- 
culation. By using equation (25), (Gr Pr) equal to 10 14 , and 
a thermal conductivity of 0.393 (Btu/(hr) (ft) (9F)) , the heat- 
transfer coefficient in the hole was found to be 15,500 (Btu/(hr) 
(sq ft)(9p)). The heat -transfer coefficient on the outside of the 
blade is usually in the neighborhood of 80 (Btu/(hr)(sq ft)(°F)). 
These values determine the temperature difference & v between the 

wall (or the hottest portion in the fluid) and the fluid entering 
the hole. By using the preceding values, it is found that for 

a large hole & v is equal to 9.1° F. Evaporation occurs when the 

wall temperature reaches or exceeds the saturation temperature of 
the water. In the holes a considerable pressure increase occurs 
with increasing distanoe from the axis of rotation. This increase 
can easily be calculated from a balanoe between the pressure and 
centrifugal forces. The saturation temperature rises with the 
pressure and it is found that an increase in radius of only 0.004 
inch is sufficient to increase the saturation temperature by 
9.1° F. Therefore, even if the liquid enters the holes in a 
saturated state, evaporation will occur only at the entrance. 

With the reduced circulation in the experimental turbine 
and the corresponding decrease in heat -transfer coefficient in the 
holes, conditions era greatly changed and evaporation may occur in 
the inndr part of the holes. 

The velocities connected with the free -convection flow are 
extremely high when the holes are large enough in diameter, which 
results in very high heat-transfer coefficients and is a definite 
advantage of the free -convection cooling as compared with forced- 
convection cooling. With the forced -convection method it is 
ordinarily impossible to obtain equally high velocities. 


FREE -CONVECTION HEAT TRANSFER IN LOOP 
Turbulent Flow 

The circulation in the holes with small diameters decreases 
because the boundary layer does not leave enough space through which 
the cold liquid may enter the holes in a direction opposite the flow 
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in the "boundary layer. It is therefore to "be expected that the 
circulation in the holes can "be materially improved when the cold 
fluid is fed into the holes from the top of the "blade (fig. 9). 

In the central part of turbine "blades, it is always possible to 
provide holes with larger diameters. If a connection is made on 
top of the "blade "between these larger holes and the small holes 
in the leading and trailing edges of the "blade, cold liquid will 
flow from the large -diameter holes to the small holes. The pro- 
posed arrangement, of cooling passages in free -convection liquid- 
cooled turbine blades in shown in figure 10. 

In "order to see how much the heat transfer in a small hole oan 
be improved in this manner, a calculation is made as follows : 

Hole (1) (fig. 9), with a very small diameter, is connected on top 
with a second hole (2) of very large diameter. The diameter of 
hole (1) may be so small that the boundary layers fill out the whole 
oross section and the liquid flows only in a downward direction. The 
velocity profile will then have practically the same shape as in 
forced -convection flow. The same thing will be true for the value 
of the friotion forces. The calculations using the friction factor 
and heat -transfer coefficient of forced flow in a tube are therefore 
applicable. First, the change of the bulk temperature t of the 
liquid along the hole may be- calculated. For this purpose, a 
heat balanoe is set up for the length dx of the tube, assuming 
the temperature t w of the tube wall to be constant along the whole 
length. Denoting the average velocity in the hole by U and the 
heat -transfer coefficient built with the local temperature difference 
ty-t-ij by Hjj yields 

CpgpU dt t = HfcrtDCt^) dx (37) 

By introducing the difference between the temperature of the 

wall and the bulk temperature of the liquid and the difference 
between the temperature of the wall and the cold liquid entering 
the hole * t w -tg, equation (37) transforms into 


d©b -4Hb 
0 b = 


CpgpDU 


( 38 ) 
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Integration with the condition that for r = 0, © b = 6 V leads 
to the expression 

gb; 

e„ . e, e ‘s° P ®u 

and for the temperature difference at the exit of the hole (x = L), 

4H-bL 

%e ' ®v « " g V ro 


The logarithmic mean temperature difference that is to "be used for 
heat -transfer calculations is 


®m “ 


6 , - e 

b,e w 


log, 


b t e 
ft 


v 


Introducing the . preceding -values gives 


e m c^gp DU 

e v " 4HjjL 



e 


gCpPDU 


(39) 


The force that generates the movement of the liquid through 
the loop is the sum of all the buoyancy forces. The pressure 
difference connected with these buoyancy forces along the length > 
dx is 


IP = gfPg-Pfc) dx 


Introducing the expansion coefficient in the same way as was pre- 
viously done and assuming the changes in speoific weight to be 
small gives 


dp - gp 6 P(t b -t 6 ) dx = gp 8 p(© v -0 b ) dx 
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The pressure difference available to overcame the friction in the 
small -diameter hole is 


(40) 


This pressure difference is used up by the friction forces. 
It is assumed that the flow is turbulent in the small -diameter 
hole and that the formula for pressure drop for turbulent forced- 
convection flow as given in reference 4 (p. 433) can, be used: 

Ap « 0.316 — ~ . p 

W -/ 4 


uf.L 

2 D 


(41) 


Ap 


dp = gp & pe v 



- e 


4Bb L 

gCppOTj 


Introducing this expression into equation 
dimensionless values gives 


0.158 



= 1 


® e D^ r D 
L 



(40) and converting into 


. 4 Wtt b 
_ e " 



(42) 


where Nu^ » HD/k. The expression on the right side of the equation 
can be simplified because it turns out, in the range of G-rashof 

numbers of 10-*- 3 to 10^ 3 and of length-diameter ratios less than 
50, that the exponent of e is a small value. The expression on 
the right side of equation (42) has the form 1 - l/z(l-e -z ). 
Transforming the e function into a series and considering only 
the first three terms gives 


1 




(43) 


The Nusselt number for forced -convection turbulent flow in a 


tube is given in reference 7 (p. 


168) as 


0.4 


NUjj = 0.023 Re-p°*®Pr 


(44) 
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Therefore, equation (42) simplifies to 

Bep 1 - 75 12.66 Hte b L _ Q.291 L 

Gr (^j 3 Ee E Pr D Eep 0 *^ 0 * 6 D 

which may be solved for the Reynolds number 


(45) 


Rep 


0.551 Gr°.515 / p\l .03 
Pp0 • 303 Xj J 


(46) 


In equation (44), the Nusselt number and the Reynolds number are 
based on the diameter of the tube and the heat -transfer coefficient 
Ejj with the mean difference between the wall and the liquid bulk 

temperature . In the previous calculations, the length x and the 
difference & v between the wall temperature and the temperature of 

the entering cold liquid was used as a basis for the heat-transfer 
coefficient and the Nusselt number. In order to make possible a 
comparison of Nusselt numbers, equation (44) is converted to this 
Nusselt number. The specific heat flow q w from the wall into the 

liquid must be the same expressed by both heat-transfer coefficients 
Therefore 


Vm * E0 v 

Inasmuch as the calculations here are for the whole tube 
length (z s* L), the Nusselt number Nu used in previous cal- 
culations is connected with the Nusselt number Su^, based on 

the diameter and the mean temperature difference, by 

*•^ 5 ^ < 47 > 


Introducing the values for the mean temperature difference © m 
from equation (39), the values for Nu^ from, equation (44), and 
simplifying according to equation (43) gives 
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Nu = 0.023 Re D 0 * 8 Pr 0 * 4 g. ML - 


0.046 


Re D °* 2 Pr 0 * 6 ^ 


(48) 


The Reynolds number may he inserted from equation (46), 


Nu = 0.0139 Gr 


0.41 pyO.154 fL\ 0 ' 176 


D/ 


1 - 


0.0522 


Qr 0.103 pyO.538 \D 


L\1 .21 

(49) 


Laminar Plow 

When the Reynolds number is smaller than the critical value of 
2100 (reference 7, p. 154), the flow in the small -diameter hole is 
laminar. Xn this case, the following formulas for the pressure drop 
and heat-transfer coefficient in laminar forced flow are used: 


Ap = 32 

Tr 

> 

Nuk ** ~ = 3.65 

* k J 


(50) 


The formula for the Nusselt number (reference 4, p. 530) is derived 
from calculations by Nusselt and is valid only far fully developed 
flow. In the inflow region, higher values occur and another more 
complicated formula has to be used. Because the simple formula (50) 
gives minimum values, it is used here. Carrying out the calcu- 
lations in the same manner as for turbulent flow results in 


32 Itejj 
Gr 


a ) 3 


i 


Re^Pr D 
14.6 L 



14.6 L' 
RejjPr D 


Nu 



Be P /L \ 3 
Gr \D / 


(51) 


(52) 
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Discussion, of Results 

The Nusselt numbers were calculated from equation (49) for 
a loop, the mna.ll hole of which has ratios L/D * 25 anil 50. The 
results are shown in figure 5. With equation (46), it was ascer- 
tained that in the range shown in figure 5 the flow is turbulent. 

The m i nim u m Reynolds number calculated was 20,000. It can. be seen 
that very high Nusselt numbers are obtained in this way. fhe Nusselt 
numb ers are even larger than in a hole closed on top with viiry large 
diameter-to-length ratio. 

The large heat-transfer coefficients obtained by means of the 
loop oan be used in turbine-blade design, a fact that will be dis- 
cussed in the last section of the report. 


ACTION OF CORIOLIS FORCES 

Up to now the action of the Coriolis forces on the flow of the 
cooling liquid In the holes has been neglected. In coder to get an 
estimate of the effect of these forces on the flow and the heat 
transfer, the magnitude of these farces is compared with the 
centrifugal forces. In the holes the relative movement of the fluid 
is normal to the circumferential velocity of the blade. In this 
case, the Coriolis acceleration Is 


2u(0 (53) 

The acceleration has the same direction as the circumferential 
velocity when the relative velocity points toward the axis of 
rotation. When the relative velooity is directed away from the 
axis of rotation, the Coriolis acceleration aots in a direction 
opposite to the circumferential velocity. The comparison with 
the centrifugal forces will be made in the cross seotion of the 
hole fixed by the mean radius r m . In a cross section of the hole, 

the relative velocities of the fluid vary between zero and the 
maximum value as previously indicated. The greatest difference 
between the Coriolis acceleration acting on different fluid particles 
In the cross section therefore has the value 


2u max co 
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^^ S ^ acoelerafcloa is com P are( ^ vith the aoceleration-oonneoted with 
the buoyanoy forces : 

r m (£pe 

The maximum difference in the temperature difference in one cross 
section is . Therefore , the maximum "buoyant acceleration due 
to the centrifugal force is 


V^ e v 

The ratio £ of tbs maxim u m Coriolis acceleration to the 
acceleration connected vith the buoyant forces is 


t - 2u max 


(54) 


The velocity may he expressed by the corresponding Reynolds 

number * Expressing, in ad d ition^ the angular velocity w by 

the Grashof nuaiber r^pe^x 3 /^ 2 , which has to "be used for con- 
vection caused by centrifugal fields, gives 


* = 2 Eeffiax (5S) 

Inserting the Reynolds number from equation (28) yields 

. 1 /x/r_\o.5 

u 0,71 ^o.m (56) 

At the mean radius r m , x is equal to L/2. The ratio x/r m has 

values near l/4. In the tests on the free -convection liquid -cooled 
turbine, the highest value of the temperature difference 0 was 
315° F. The film temperature for evaluating property values of the 
cooling water was 200° F. When this film temperature is used, the 
Prandtl number is 1.74 and the expansion coefficient 0.00043/° F. 
Inserting these values into equation (56) gives the value 0.84 for 
the ratio of the Coriolis acceleration to the acceleration connected 
with the buoyancy forces in the cross section of the hole . The 
Coriolis forces are therefore quite appreciable. 
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The Coriolis forces set up a secondary movement of the type 
shown, in figure 11. The secondary flow is similar to tbs one in 
a hent tube. It increases the boundary-layer thickness on the side 
of the hole in the direction of the rotation and decreases its 
thickness on the opposite side. 

To some extent it 1 b possible to study the influence of this 
secondary movement on heat transfer by means of static tests pre- 
viously described. By inclining the static test section, the free- 
convection-flow setup would be exactly similar to the flow in the 
rotating hole in the turbine blade if the velocity and temperature 
profiles would be similar. In reality, the profiles differ In 
shape; however, the comparison may be used as a first estimate of 
the Coriolis forces. The ratio of Coriolis to centrifugal accele- 
ration of 1.00 corresponds to a static tube Inclined at 45° from 
the vertical. Unfortunately, no test results are known for gravita- 
tional free -convection flow in an Inclined heated tube. Measurements 
were made, however, on gravitational free -conveot ion flow in an air 
layer between two parallel plates, vertical and inclined, and in an 
air layer between two coaxial cylinders. In both cases It was found 
that the heat -transfer coefficient changed very little when the 
test setup was Inclined from the vertical position to an angle of 
45°; the change was about 20 percent (reference 4) . This infor- 
mation gives some indication that even when the ratio of the 
Coriolis to the centrifugal forces Is considerable, the ratio does 
not affect the heat transfer to a large degree. 


SUMMARIZING REMARKS AND RECOMMENDATIONS 

The method of free -conveotion cooling of rotating turbine 
blades utilizes the incense .currents generated in the liquid, which 
fills the hollow spaces within the blades, by the centrifugal farces 
as soon as temperature differences are present In the liquid. In 
the current designs using this method, holes are provided within 
the blades that are closed on top and connected with some liquid 
supply at the blade root. Such a blade was shown and the free- 
convection currents set up in one hole of the blade were illustrated. 
The heated fluid within the boundary layer along the wall of the hole 
moveB toward the blade root and. is replaced by cold liquid flowing 
toward the top in the central part of the hole. When the diameter 
of the hole is very fwn*n as compared with the length, the two liquid 
streams flowing in opposite directions hinder each other and the 
circulation and oooling effect deorease. A measure of this limiting 
diameter-to-length ratio was obtained by calculating the boundary- 
layer thickness of the free -conveotion flow in a hole with a very 
large diameter-to-length ratio. 
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Inasmuch as the flov in the boundary layer is turbulent under 
the conditions present in a rotating turbine blade, the approximate 
method introduced by von Karmdn vas used for the calculations. In 
this method the shapes of the velocity and temperature profiles 
have to be assumed. Unfortunately, no information is available on 
the shape of either profile in turbulent free-oonvection flow. The 
calculations were therefore carried through with three different 
shapes for the velocity profile. The shape of the temperature pro- 
file was assumed to be identical to the one in turbulent forced - 
convection flow. Also, the law for the shearing stress on the wall 
in forced flow was used. The method proved quite insensitive to a 
change in the shape of the profiles. The values of the boundary- 
layer thickness and of the maximum velocity for the flow in the 
range covered by the results of the calculations with the different 
assumptions varied approximate ly 20 percent. A more accurate 
answer to the problem can be obtained only when more experimental 
information is available on the shapes of the profiles. 

The calculations and their comparison with test results on a 
free -convention cooled turbine indicated that in blades of 2.5-lnch 
length the circulation and the cooling effect are decreased to a 
considerable degree in holes with diameters less than 0.146 inch. 

The circulation and heat removal of a hole with a small diameter -to- 
length ratio can, however, be restored to the value of one with a 
large diameter-to-length ratio when the hole with the small, diameter- 
to-length ratio is connected on top of the blade with neighboring 
large -diameter-to-length-ratio holes. A second set of calculations 
proved this fact. 

Prom the calculations, the following recommendations for the 
arrangement of the oooling holes in turbine blades can be made. The 
blade should be thick in the central part so that holes with large 
diameters can be arranged there. These holes should be oonneoted 
on top of the blade with the small -diameter holes. Such an arrange- 
ment of oooling holes in a turbine blade was shown. For especially 
effective oooling, the cross seotion of the holes should be as large 
as possible. It must be kept in mind, of course, that very high 
pressures are present inside the holes and that the walls must with- 
stand these pressures. The stresses caused by the internal pressures 
in the holes must be added to the stresses caused by the centrifugal 
forces in the solid material of the blades. 
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In tiie boundary-layer calculation, the effect of the Coriolis 
forces on the movement of the liquid in the rotating holes was 
neglected. Experiments on gravitational free -oonvect ion flow on 
inclined objects indicated that this effeot is of little importance 
on the heat transfer. 


Lewis Flight Propulsion Laboratory, 

national Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

C constant 

C s constant for variation of heat-transfer coefficient with x 1 / 5 

C constant for variation of maximum velocity in boundary layer 

with x m 

C g constant for variation of boundary-layer thickness with x n 

Cp specifio heat at constant pressure, Btu/(lb)(°P) 

D diameter, ft 

Gr Grashof number, 

u2 

g acceleration due to gravity, ft/sec^ 

H heat -transfer coefficient, Btu/(sq ft)(sec)(°F) 

k heat conductivity, Btu/(ft) (sec) (°F) 

L length of hole, ft 

l length, ft, (fig. 2) 

m exponent 

Nu Nusselt number, — 

’ k 

n exponent 

U gO-JX 

Pr Prandtl number, — = 

7 a k 

p pressure, Ib/sq ft 

q specific heat flow, Btu/(hr)(sq ft) 

Eejj Eeynolds number based on diameter D, SL 
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® e max 

r 


t 

tr 


u 

U 1 


V 

X 

y 

a 

P 

5 

8 * 

i 


e 

■ 

H- 

V 

P 

T 

CO 


Eeynolds number 'based on maxi, mum velocity -2|p — 

radius, ft , 

temperature, °F 
velocity, ft/seo 

velocity component In x -direction, ft/seo 

velocity outside boundary layer of comparable forced- 
convection flow, ffc/seo 

specific volume, cu ft /lb 

coordinate (distance from stenting point of boundary layer), ft 

coordinate (distance from weill), ft 

thermal diffusitivity, set ft /sec 

expansion coefficient, l/°F 

boundary-layer thickness, ft 

displacement thickness of boundary layer, ft 

ratio of maximum Coriolis acceleration to TnaTlnmTa oentrifugal 
acce" 1 '“ration 

temperature difference, °P 

temperature difference between wall and fluid outside of 
boundary-layer, °F 

dynamic visoosity, (lb) (sec)/sq, ft 

kinematic viscosity, sq. ft/sec 

dynamic density, (lb) (sec^J/ft^ 

shearing stress, lb/sq ft 

angular velocity, l/sec 


Subscripts : 


b 


bulk value (mixed mean) 
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e exit condition 

m mean value 

max maximum value 

v on vail 

5 on outer border (and outside) boundary layer 

I values calculated with first assumption for velocity profile 

II values calculated with second assumption for velocity profile 
ITT values calculated with third assumption for velocity profile 



1307 


NACA EM E5QD25 


35 


APPENDIX B 


SOLUTION OF EQUATION (12) 
The differential equation, to he salved is 


• ^“max „ sP 0 v %aT 

— = 0.254 - — - — 

QJC U 3HELX * 



"max 2 * 


1/4 


^ 5/4 


( 12 ) 


It is advantageous to transform this equation by the use of non- 
dimensional values. The Grashof number Gr is used to replace 

the independent variable x and the Eeynolds number Be ^^ * Um ^ xX 
for the dependent variable Differentiating the Eeynolds 

number gives 


a ( Ee max) ^ ^ du max . ^max 
dx v dx v 


(57) 


Introducing the differential quotient dUg^/dx from this expression, 
and the Grashof number into equation (12) yields 


d(Ee ) _ Q.Q 852 _ 0.00194 Be^ 2 Pr 5 / 4 
a ( Gr > " Ee max Gr Nu 5 A 

The Nusselt number is replaoed by use of equation (10), which gives 

4®*^) 0.0652 °- 0250 E W 2 SrS/6 

d(Gr) = R emax Gr 17/12 {SB) 

This equation was graphically solved (fig. 12). Inasmuch as the 
range of Grashof numbers concerned is fairly wide, it is advantageous 
to obtain the graphical solution in a logarithmic diagram in which 
the Eeynolds number is plotted over the Grashof nu m b er on a log- 
arithmic scale. In order to determine the curves that present the 
solution of the differential equation in the diagram, the slope for 
any pair of values Be and Gr is needed. If logg Be ma T = and 

logg Gr = |, then the . needed slope is 
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dJj/ B Gr 

d| Eemax 4 (G**) 


For any pair of -values Rejnax 6114 Gr > bhe differential quotient 


dCRe^j) 

4(Gr) 


can be determined from equation (58) . 


Equation (59) 


gives the slope in the logarlthmlo plot. Figure 12 shows the slopes 
determined in this way for two Frandtl numbers . The needed solution 
of the differential equation (58) must fulfill the condition that 
for Gr = 0, R e mai “ 0» This condition follows from the assumption 

that the turbulent boundary layer starts with a thickness zero at x 
equal to zero. The validity of this assumption was previously dis- 
cussed. In starting out with high Grashof numbers and drawing curves 
that have the described slopes for all Grashof numbers, it is found 
that the curve that fulfills the preceding condition can be found 
with very good accuracy-. All ourves that start out with a Reynolds 
number that is too great turn away in an upward direction as the 
Grashof number is decreased. Curves that start cut at high Grashof 
numbers with a low Reynolds number turn away in a downward direction. 
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Figure 6. - Calculated displacement thickness 6* of turbulent f ree-convect i on 
boundary layer plotted against product of Grashof number and Prandtl number. 
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Figure 7. - Calculated displacement thickness 6* of turbulent f ree-convect i on 
boundary layer plotted against Prandtl number, H proportional to x^ 5 ; 

(Gr Pr) - |0 12 . 
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Figure 8. - Reynolds number Re raax characterizing maximum velocity In boundary layer and displacement 
thickness fi plotted against product of Grashof number and Prandtl number calculated with two 
assumptions for heat-transfer coefficient. Prandtl number = 1.0. 
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Figure 12. *■ Graphical solution of differential equation 
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